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1. Introduction 



Nucleon sigma term o^n is given by a scalar form factor of nucleon at zero recoil. While up 
and down quarks contribute to o n N both as valence and sea quarks, strange quark appears only as 
a sea quark contribution. As a measure of the strange quark content of the nucleon, the y or fj s 
parameters are commonly introduced. These parameters are defined as 

/mi- j fi»r\ 2{N\ss\N) , m s (N\ss\N) 

a JzN =m ud (N\uu + dd\N), y= ]_ - ' f T = ' (1.1) 

{N\uu+dd\N) M N 

The y parameter plays an important role to determine the detection rate of possible neutralino dark 
matter in the supersymmetric extension of the Standard Model [JT], ||]. 

Using lattice QCD, one can calculate the nucleon sigma term directly. Furthermore, it is pos- 
sible to determine the valence and sea quark contributions separately. Previous lattice studies of the 
nucleon sigma term and strange quark content have been done within the quenched approximation 
[g, f§], or a two-flavor QCD calculation [g, g|. Recent studies use 2+1 flavor QCD tffy. 

There was an apparent puzzle in these results: the strange quark content is unnaturally large 
compared to the up and down contributions that contain the connected diagrams too. Concerning 
this problem, it was pointed out that using the Wilson-type fermions, the sea quark mass depen- 
dence of the additive mass renormalization and lattice spacing can give rise to a significant lattice 
artifacts in the sea quark content [§]. Our previous study from a two-flavor QCD simulation using 
the overlap fermion removed this problem by explicitly maintaining exact chiral symmetry on the 
lattice for the first time. It was revealed that the sea quark gives only a small contribution to the 
nucleon sigma term. 

In this study, we analyze the data of the nucleon mass obtained from a 2+1-flavor QCD simula- 
tion employing the overlap fermion We employ the chiral perturbation theory (ChPT) analysis 
for the nucleon mass data and calculate the sigma term from its up and down quark mass depen- 
dence. To calculate the strange quark mass dependence, we use the reweighting technique for 
the strange quark mass. Our study with exact chiral symmetry avoids the contamination due to a 
significant lattice artifact so that it provides a reliable calculation of the strange quark content. 



2. Lattice simulations 



We make an analysis of the nucleon mass obtained on 2+1-flavor QCD configurations gener- 
ated with dynamical overlap fermions. Our simulations are performed on a 16 3 x 48 lattice with a 
trivial topological sector Q = 0. For the gluon part, the Iwasaki action is used at /? = 2.30 together 
with unphysical heavy Wilson fermions and associated twisted-mass ghosts [[R]] introduced to 
suppress unphysical near-zero modes of Hw(— mo). The lattice spacing a = 0.108 fm (a -1 = 1.83 
GeV) is determined through the Sommer scale ro = 0.49 fm of the static quark potential. We 
take two different strange sea quark masses am s = 0.080 and 0. 100 with five values of degener- 
ate up and down quark masses (for am s = 0.080, am sea = 0.015,0.025,0.035,0.050,0.080 and for 
am s = 0.100, am sea = 0.015,0.025,0.035,0.050,0.100). We accumulate about 2,500 trajectories 
for each combination of sea quark masses; the calculation of the nucleon mass is carried out at ev- 
ery five trajectories, thus we have 500 samples for each sea quark masses. In order to improve the 
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Figure 1: (Left) The nucleon mass data for the partially quenched data points. The empty and full symbols 
show the fit results for am s = 0.080 and am s = 0.100, respectively. (Right) Chiral fit of the lattice raw 
data(full symbols) and the finite size corrected data(empty symbols). The solid and dashed curve represent 
the fits of raw and corrected data. For a reference, we also show the experimental values of the nucleon mass 
(star). 

statistical accuracy, we use the low-mode averaging technique. The two-point function made of the 
low-lying modes of the Dirac operator is averaged over different times slices with 80 chiral pairs 
of the low modes. For the quark propagator, we take a smeared-local two-point correlator and fit 
the data with a single exponential function after averaging over forward and backward propagating 
states in time. The statistical error is estimated using the starndard jackknife method with a bin 
size of 10 samples. We take the valence quark masses am m i = 0.015, 0.025, 0.035, 0.050, 0.060, 
0.070, 0.080, 0.090 and 0. 100. All the plots of the nucleon mass spectrums are shown in the left 
panel of Fig. [j] 

The matrix element defining the nucleon sigma term can be related to the quark mass de- 
pendence of the nucleon mass using the Feynman-Hellman theorem. For the strange sea quark, the 
matrix element is represented by the following relation as = (A^l^ |A^)disc - The subscript "disc" 
on the expectation value indicates that only the disconnected quark line contractions are evaluated. 
In the present study we exploit this indirect method to extract the matrix elements corresponding 
to the nucleon sigma term. 

Another possible method is to directly calculate the matrix element from three-point functions 
with an insertion of the scalar density operator (uu + dd){x) or ss[x). We note that the indirect 
method gives mathematically equivalent results to those by the direct method including lattice arti- 
facts. Numerical differences could arise only in the statistical error and the systematic uncertainties 
due to the fit ansatz. 

3. Analysis of the nucleon mass using chiral perturbation theory 

We carry out a fit with the SU(2) NLO ChPT for the baryon, 

3„2 

M N = M Q -4c l ml-^^ml + e' l (n)m 4 Jt , (3.1) 

where Mq is the nucleon mass in the chiral limit and f n is the pion decay constant fixed at its 
physical value 92.4 MeV and the constant gA describes the nucleon axial-vector coupling. We fix 
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without FVCs with FVCs 

4 points 5 points 5 points 

am s = 0.080 50.0(4.5) 48.2(3.0) 51.1(3.0) 

am s = 0.100 51.3(4.1) 45.5(3.0) 48.2(3.0) 

Table 1: Nucleon sigma term g k n [MeV] with and without the finite volume corrections (FVCs). The error 
is statistical only. 



gA to its experimental value g& = 1.267. We use the pion mass obtained from the calculation of the 
meson spectra. Since our lattice volume is not large enough, the finite size effects must be taken 
into account especially in the data at lightest quark mass. To estimate possible systematic error, we 
attempt the fit with and without lightest sea quark mass (4 points and 5 points, respectively). We 
also correct the finite volume effect based on the NLO ChPT formula. The nucleon mass difference 
in the finite volume with box size L and infinite volume 8M^{L) is calculated within the ChPT. We 



use the asymptotic formula Qll| ] to evaluate the magnitude of the finite volume effect as 



5Mn ^ = f St + ^r¥r^ 1 7 ex p(-^)> ^ 

\ lOTT/w (27T)lAL2 / L 

where A is the A-nucleon mass splitting, and the coupling g^N is decuplet-octet-axial coupling. We 
use the phenomenological estimates g&N = 1-5 and A = 300 MeV. 



The chiral fit is made for the corrected five data points using the formula (p.l|). The results are 
shown in the right panel of Fig. [T[ After correcting the finite volume effect, there are 6% decrease in 
Mq and 5% increase in the magnitude of the slope \c\ |. Finally we obtain the nucleon sigma terms 
by estimating the pion mass dependence of the nucleon mass at the physical pion mass. The results 
are shown in Table [jj We observe that both the fit results at am s = 0.080 and am s = 0. 100 without 
lightest points are consistent each other. These results are also consistent with our previous results 
for two-flavor simulations. We quote our best value of the sigma term as a fit result of 4 points 
fit (am s = 0.080, FVCs not included). We obtain o k n = 50.0(4.5). The finite volume correction 
of the sigma term amounts to about 6%. We note that a sizable chiral extrapolation error should 
also be added after comparing NLO and NNLO in ChPT as found in our previous analysis with 
two-flavor QCD. 

4. Reweighting the strange quark mass 

In order to study the strange quark mass dependence of the nucleon, we utilize the reweighting 



technique [|12j, |13Q . Using this method we can calculate any physical observables at a slightly 
different strange quark mass (m' s ) from the configurations (U) generated at a certain value of strange 
quark mass (m s ) as 

(0(U,D(m' s ))w(m>,m s )) 
\v) m i = — — ■ (4.1) 

{w(m' s ,m s )) 

where D is overlap Dirac operator and we introduce the reweighting factor w(m' s ,m s ) defined as 



w(ra' mJ] = det 



£K) 
OK) 



. We calculate the reweighting factor by decomposing it into the low and 
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trajectory 

Figure 2: The fluctuations of the reweighting factor. The plot shows the noise number dependence of the 
fluctuations at am va i = am sea = 0.050 with strange quark masses am' s = 0.075, am s = 0.080. 

high-mode contributions. A former contribution is obtained exactly by the 80 pairs of low-lying 
eigenvalues of the Dirac operator. For high-mode, we use the improved stochastic estimator 

w high (m>,) = ^l^-K^o-D^, (4.2) 

a = D{m s f{D{m' s )- i yD{m' s )- l D{m s ), 

where B, r are the ensembles of pseudo-fermion fields (N r : number of noise vector), and each are 
generated with Gaussian probability. Here we also introduce the projection operators P and P to 
project the eigenspaces of £ into low and high-modes, respectively. 

Fig. ^| shows the reweighting factor with am' s = 0.075 and am s = 0.080 for each gauge con- 
figuration. It is evaluated with different number of noise vector N r =0, 5, 10 and 50. We find that 
the statistical error of the reweighting factor is very small compared to that from the fluctuation 
of the gauge configurations. We also find that the contribution from high-mode is almost indepen- 
dent of the gauge configuration and the configuration dependence is almost completely dictated by 
the low-mode contribution. These properties are universal for entire region of up and down quark 
masses. 

In order to estimate the strange quark mass dependence of the nucleon mass, we evaluate the 
nucleon masses at 10 different strange quark mass for each configurations. That covers the strange 
quark mass range am' s = 0.0775-0.0825 for am s = 0.080 and am' s = 0.0975-0. 1025 for am s = 0.100. 
To compute the high-mode contributions for several m' s at once, we use the multishift solver for the 
Dirac operator. We take five random noises for the noisy estimator. Thus the nucleon mass for each 
different strange quark masses is calculated. We find small strange quark mass dependence of the 
nucleon mass and its statistical error does not grow for entire region of the strange quark mass we 
take. 

5. Extraction of the strange quark content 

To extract the strange quark content (AT|.f.y|Af), we fit the nucleon mass by a linear function of 
strange quark mass. Then the coefficient of the linear term corresponds to the strange quark content 
of the nucleon. We use 1 1 data points for the fit including the unreweighted point at each sea and 
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Figure 3: The chiral extrapolation fit for the strange quark content for m s = 0.080 (black) and m s — 
0.100 (red). The plot shows only the data at unitary point. The linear line shows the fit functions for 
unitary point. For a reference, the empty symbols show the data not used in the fit. 



am sea = am va i 


0.015 


0.025 


0.035 


0.050 


0.080 


0.100 


am s = 0.080 
am s = 0.100 


0.25(28) 
0.04(20) 


0.30(10) 
0.15(15) 


0.31(15) 
0.16(13) 


0.42(11) 
0.22(10) 


0.01(7) 


0.41(10) 



Table 2: The strange content of the nucleon at the unitary points. 




valence quark masses. We obtain the non-zero signals of the strange content at each data point 
with 50-100% statistical error. The results for the unitary points are shown in Table ||. To estimate 
the physical values of the strange content, we need to carry out a chiral extrapolation. In order to 
improve the accuracy, we use partially quenched data point. We take the fit ranges as am sea = 0.025- 
0.050 and am va i = 0.025-0.050 to avoid large finite volume corrections and possible higher order 
corrections. Its chiral behavior is described by a derivative of the nucleon mass formula in SU(3) 
ChPT in terms of the strange quark mass. We note that the non-analytic behavior of the light up 
and down quark masses does not appear up to NNLO corrections. Therefore we safely use the 
polynomial function of the up and down quark masses. We use a simple fit function 

dMN , ,< n 

— — =a\ + a 2 m sea + a 3 m va i. (5. 1) 

am s 



The fit results are shown in Fig. |5[ The data are reasonably fitted by the formula (5J_). Although 



our results at chiral limit have about 100% errors, it shows that the physical result of the strange 
quark content is very tiny. We estimate the renormalization group invariant quantity fj s as 

= f 0.020(34) (am, = 0.080) 
lTs [ 0.019(53) (am, = 0.100) ' 

where the error is statistical. There may be some systematic uncertainties coming from the finite 
volume corrections and chiral extrapolations. Both results are consistent each other, we find the 
strange quark mass dependence of this quantity is very small. Using two results of the sigma 
term obtained in Sec. |3] and fj s (=0.020), we obtain y ~ 0.03. These results are comparable to 
our previous two-flavor analysis [Qj. We also conclude that the upper bound on fj s is given by 
|/ r J <0.08 at 1 a level. 
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6. Summary 

We calculated the nucleon mass and sigma term in 2+1 -flavor QCD simulation on the lattice 
with exact chiral symmetry. To evaluate the strange quark content of the nucleon, we used the 
reweighting technique for the strange quark mass, worked very well for studying the strange quark 
mass dependence. Although our results have large statistical uncertainty, we conclude that the 
strange quark content of nucleon is tiny compared to its up and down quark contributions. 

The main numerical calculations were performed on IBM System Blue Bene Solution at 
High Energy Accelator Organization (KEK) under support of its Large Scale Simulation Pro- 
gram (No.09-05). We also used NEC SX-8 at Yukawa Institute for Theoretical Physics (YITP), Ky- 
oto University and at Research Center for Nuclear Physics (RCNP), Osaka University. This work 
is supported in part by the Grant-in-Aid of the Ministry of Education (Nos. 21-897, 19540286, 
19740160, 20105001, 20105002, 20105003, 20105005, 20340047, 21674002, 21684013). The 
work of HF was supported by the Global COE program of Nagoya University "QFPU" from JSPS 
andMEXT of Japan. 
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